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^The  major  objectives  of  this  project  are  to  develop  computational 
algorithms,  both  sequential  and  parallel,  for  several  important  linear  algebra 
problems  that  arise  in  the  design  and  analysis  of  linear  control  systems 
governed  by  the  systems  of  ordinary  differential  equations.  include 

.  (i)  controllability  problems,  (ii)  stability  and  inertia  problems,  tlTl)  pole 
assignment  problems,‘^;^(4v)^  matrix  equations  problems,  (vX  relative  primeness 
of  polynomials  and  matrices  and  the  Cauchy-index  problems  of  rational 
functions  etc.  Besides,  a  part  of  the  project  is  devoted  to  the 
theoretical  study  of  the  parallel  arithmetic  complexity  of  these  problems*, 
that  is,  how  fast  these  problems  can  be  solved  in  parallel  assuming  that 
sufficiently  many  processors  are  available. 

Though  some  numerically  viable  sequential  algorithms  have  been  designed 
for  some  of  these  problems  in  recent  years,  parallel  algorithms  and 
algorithms  for  large  scale  problems  are  virtually  non-existent.  ^ 

I  am  very  pleased  to  report  that  we  have  taken  a  "LEAD”  in  this  area. 
Since  my^INTERIM  REPORT"  last  October,  the  following  additional 
accomplishments  have  been  made: 

(A)  Linear-time  steps  algorithms  requiring  only  0(n)  processors  for 
implementations  have  been  developed  for  controllability,  observability, 
stability,  inertia  and  Cauchy  index  problems  for  polynomials.  These  are 
sharp  improvements  over  our  0(n  log2  c)  steps  -  0(n  )  processors 
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algorithms  reported  last  year.  Especially,  savings  in  number  of 


processors  is  significant.  The  results  have  been  presented  at  the  SIAH 
Conference  on  Parallel  processing  for  Scientific  Computations,  Norfolk, 
Virgina,  November  1985.  The  paper  based  on  these  results  is  presently 
being  written  up. 

(B)  Fast  sequential  and  parallel  algorithms  have  been  developed  for  (i) 
determining  relative  primeness  and  the  number  of  common  eigenvalues  between 
two  given  matrices  and  (ii)  finding  the  inertia  and  stability  of  a  non- 
hermitian  matrix.  Sequentially,  these  algorithms  are  more  efficient  than 
the  existing  procedures  and  in  parallel,  they  can  be  implemented  in  almost 
linear  time  steps  using  at  most  0(n  )  processors.  A  paper  based  on  these 
results  has  just  appeared  in  a  special  issue  of  CONTEMPORARY  MATHEMATICS 
(Volume  47)  on  '^Linear  Algebra  and  its  Role  in  Systems  Theory**  published  by 
the  American  Math  Society. 

(C)  Jointly  with  A.  Sameh  of  University  of  Illinois,  a  parallel 
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algorithm  for  the  single-input  pole  assignment,  seemed  to  be  suitable  for  a 
shared  memory  multiprocessor  like  Cray  X  -  HP/48,  has  been  developed. 

The  algorithm  is  presently  being  implemented  on  **Alltant**  system  at 
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University  of  Illinois. 

(D)  It  has  been  shown  that  the  parallel  arithmetic  complexity  of  the 
OR  factorization  of  a  matrix  A  is  upper  bound  by  Odog^  n)  steps.  As 
an  application  of  this  result,  a  parallel  algorithm  requiring  only  0(log2  n) 
steps  and  0(n)  processors  for  computing  the  eigenvalues  of  a  symmetric 
tridiagonal  matrix  has  been  developed.  The  algorithm  has  been  implemented 
on  HEP  at  Argonne  National  Laboratory  in  an  M.Sc.  thesis  at  Northern 
Illinois  University  by  Ava  Chon  and,  the  results  were  compared  with  parallel 
Q-R  algorithm  of  Sameh  and  Kock  and  the  method  of  mult isect toning.  The 
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results  will  soon  appear  in  a  paper  in  INTERNATIONAL  JOURNAL  OF  COMPUTER 
MATHEMATICS. 

(E)  A  new  result  on  the  existence,  uniqueness  and  nonsingualirity  of 
solutions  of  a  special  type  of  Sylvester's  equation:  AX  +  XB  =  C,  has  been 
obtained.  As  applications,  direct  methods  have  been  developed  for 
constructing  symmetr izers,  ccmputing  the  characteristic  polynomial  of  a 
matrix  and  finding  the  numbers  of  common  eigenvalues  between  two  matrices. 
The  results  were  presented  at  the  SIAM  Conference  on  Applied  Linear  Algera, 
Raleigh,  NC,  1985,  and  the  paper  has  been  accepted  for  publication  in 
LINEAR  ALGEBRA  AND  ITS  APPLICATIONS. 

(F)  Two  new  algorithsm  have  been  developed  to  compute  the 
characteristic  polynomial  of  a  Hessenberg  matrix,  they  are  more  efficient 
than  the  existing  ones:  one  by  the  P.I.  and  the  other  due  to  Wilkinson. 
Also  it  has  been  shown  that  the  parallel  arithmetic  complexity  of  this 
problem  is  upper  bounded  by  OClog^  n)  parallel  steps.  The  paper  has 
appeared  in  J.  INDUSTRIAL  MATH.  SOCIETY. 

(G)  A  new  efficient  (0(n  log2  n)  steps)  parallel  algorithm  has  been 
developed  to  compute  the  zeros  of  a  polynomial.  The  paper  is  being  written 
up. 

(H)  An  efficient  highly  parallel  algorithm  has  been  designed  tosolve 
the  matrix  equation  problem:  AX  +  XB  =  C,  in  the  case  A  is  a  symmetric 

s 

matrix.  Attempt  is  being  made  to  extend  the  algorithm  to  the  nonsymmetric 
case.  The  preliminary  results  were  presented  at  the  SIAM  Conference  on 
Applied  Linear  Algelbra,  Raleigh,  NC,  April  1985. 

(I)  By  means  of  constructive  procedures  it  has  been  shown  that  most  of 

y 

the  control  problems  under  discussion  can  be  solved  in  as  fast  as  Odog^  n) 
parallel  steps  using  sufficiently  many  processors.  These  results  are 


of  theoretical  interests  only  at  present.  The  results  will  appear  in  the 
paper  "On  Parallel  Arithmetic  Complexities  of  Some  Linear  Algebra  Problems 
in  Control  Theory*'. 

A  brief  overview  of  these  results  were  presented  by  the  principal 
investigator  in  an  INVITED  session  on  "Linear  Algebra  and  Systems  Theory 
II''  at  the  International  Symposinm  on  Mathematical  Theory  of  Networks  and 
Systems,  Stockholm,  Sweden,  Tune  1985.  An  INVITED  paper  entitled  ' 
Theroetiical  and  Computational  Aspects  of  Some  Linear  Algebra  Problems  in 
Control  Theory"  based  on  this  talk  will  soon  appear  in  the  PROCEEDINGS  OF 
THE  CONFERENCE  published  by  North-Holland:  C.  Byrnes  and  A.  Lindquist 
editors.  An  one-hour  KEY-NOTE  address  was  also  given  by  the  principal 
investigator  at  the  International  Conference  on  Linear  Algebra  and  its 
Application  in  Coimbra,  Portugal,  March,  1985,  based  on  some  of  these 
results.  Besides  these  conference  talks,  several  INVITED  colloquium  talks 
have  been  delivered  by  the  principal  investigator  in  and  outside  the  Untied 
States.  They  include  colloquium  talks  at:  1)  Yale  University,  New  Haven, 
Connecticut,  2)  University  of  Illinois  at  Urbana-Champaign,  3)  Kent  State 
University,  Kent,  Ohio,  4)  Indian  Institute  of  Technology,  Bombay,  India, 

5)  Indian  Statistical  Institute,  Calcutta,  India  and  6)  Tata  Institute  of 
Fundamental  Research,  Bombay,  India,  etc. 

The  principal  investigator  will  also  present  a  half-hour  talk  at  the 
IEEE  CDC  Conference  in  Florida,  December  1985  and  has  accepted  invitations 
to  give  INVITED  talks  at  the  (i)  International  Conference  on  Linear  Algebra 
and  Analysis,  Czechoslovakia,  March,  1986,  and  (ii)  American  Control 
Conference,  Seattle,  Washington,  June  1986. 

In  designing  the  parallel  algorithms  reported  in  this  paper,  a  special 
attention  has  been  given  to  the  use  of  the  existing  efficient  parallel 


algorithms  for  various  types  of  matrix  computations.  As  a  result,  these 
algorithms  have  potentials  for  immediate  implementations  on  some  of  the 
existing  and  presently^bnilt  parallel  machines.  Our  present  attention  is 
being  focused  on  implementat ing  some  of  the  algorithms  on  DAP,  Alliant 
System  at  University  of  Illinois,  Urbana-Champaign,  and  other  available 
machines,  improve  and  redesign  the  already-developed  algorithms  for  the 
existing  machine.s  and  develop  algorithms  for  Large-Scale  Control  problems. 
In  the  area  of  Large-Scale  computations  with  control  problem,  a  lot  remains 
to  be  done  in  yearrs  ahead  -  the  field  is  wide-open. 
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Linear-time  step  algorithms  requiring  only  0(n)  processors  for  implementation  have  been 
developed  for  controllability,  observability,  stability,  inertia  and  Cauchy  index  problems 
for  polynomials.  These  results  were  presented  at  the  SIAM  Conference  of  Parallel  Processing 
for  Scientific  Computations  in  November  1985.  Also,  a  paper  on  fast  sequential  and 
parallel  algorithms  for  determining  relative  promeness  and  the  number  of  common  eigenvalues 
between  two  given  matrices  and  finding  the  inertia  and  stability  of  a  nonhermitian 
matrix  has  appeared  in  a  special  issue  of  Contemporary  Mathematics  on  "Linear  algebra  and 
its  role  in  systems  theory". 
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